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Abstract 

Advances in first-principles computational approaches have, over the past fifteen 
years, made possible the investigation of physical properties of ferroelectric sys- 
tems. In particular, such approaches have led to a microscopic understanding of 
the occurrence of ferroelectricity in perovskite oxides at ambient pressure. In this 
paper, we report ab-initio studies on the effect of hydrostatic pressure on the fer- 
roelectricity in perovskites and related materials. We found that, unlike commonly 
believed, these materials exhibit ferroelectricity at high enough pressure. We also 
explained in details the (unusual) nature of this ferroelectricity. 
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I. INTRODUCTION. 



ABO3 perovskites form one of the most important classes of materials because they can 
exhibit a broad range of properties, e.g. superconductivity, magnetism, ferroelectricity, 
piezoelectricity, dielectricity and multiferroism. Such properties can be varied - and thus 
optimized to generate various devices with great performance and diverse functionalities 
- thanks to many factors. Examples of such factors are external magnetic and electrical 
fields, atomic substitution, chemical ordering and pressure. 



In a classic paper ((Samara et all Il975l ). Samara at al. were able to explain the effect of 
hydrostatic pressure on displacive phase transitions associated with soft zone-center trans- 
verse optic (TO) and zone-boundary phonons. Among the most striking pressure effects 
is the decrea se of the transition t emperature with pressure and ultimate vanishing of fer- 



roelectricity ((Samara et al 



1975). As it was shown, the vanishing of fer roelectricity with 



pressure can be readily understood from the well-known soft-mode theory (jOinzburg 



Mill) 



where the soft mode frequency depends on the short-range and Coulomb interactions. In 
this case, the vanishing of ferroelectricity is associated with a much more rapid increase of 
the short-range interactions than the long-range interactions with increasing pressure. As a 
result, the harmonic soft-mode frequency becomes less negative with increasing pressure. It 
is interesting to note that for the zone boundary modes the roles of the short-range and long- 
range forces are reversed, which leads to an enhancement of these modes when increasing 



pressure. Despite its simplicity, this mode l could successful 



of many compounds ((Samara et al 



1975 



Sani et al. 



y explain high pressure behavior 



2004|). [Note that "long-range" inter- 



actions denote interactions that persist even at distances where the electronic wavefunctions 
no longer overlap. Such long-range interactions have a distance dependency that is inversely 
proportional to some power of the interionic distance. "Short-range" interactions only occur 
when the electronic wavefunctions overlap. They typically have a distance dependency that 
follows an exponential form at large distances, and are inversely proportional to some power 
of the interionic distance at smaller distances. 

The microscopic understanding of the occurrence of ferroelectricity in perovskite oxides 
has been substantially improve d by first-principles electronic band-structure t echniques. In 



particular, Cohen & Krakauer ((Cohen 



1992; 



Cohen and Krakauer . 



1990 



1992) showed that 
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a delicate balance exists between Coulomb interactions that favor ferroelectric distortions 
and the short-range repulsions that favor the undistorted high-symmetry structure. It it has 
been established that the balance between different competing interactions can be tipped by 
relatively small effects that are related to covalency, namely hybridization between p orbitals 
of the oxygen atoms and the d orbitals of the B atoms. In fact, such covalency is essential 
for stabilizing the ferroelectric distortion or, equivalently, triggers the soft-mode instability. 
The reason for this stabilization is that the hybridization weak ens the short-range repu lsions. 
This fact has been elegantly demonstrated by Posternak et al. ([Posternak et a/.Ul995T ) - who 
showed that the ferroelectric instability disappears in KNbOs when the interaction between 
the O 2p and Nb 4d states is artificially suppressed. Since then, similar conclusions have 
been drawn for other ferroelectric perovskites. Moreover, an additional Pb-0 hybridization 



is the ke y factor to explain why PbT i O.s has deeper 



BaTiOs ((Cohen and Krakauer 



1992: 



Kuroiwa et al. 



erroe 



ectric double-well potentials than 



20011 ). The important role of covalent 



bonding is also clearly reflected in the anoma lously large values of the Born dynamical 



charges (jChosez et al 



1998; 



Zhong et al 



19941 ) 



In this paper, we use first principles calculations to investigate many insulating per- 
ovskites and related materials under hydrostatic pressure, and show that the conclusion 
about the role of the different interactions menti oned above is indee d well justified at low 



pressure but is no longer valid at high pressure ([Kornev et al 



20051 ). More precisely, our 



study reveals (1) the absence of a cubic paraelectric phase in favor of a polar phase under 
high-enough hydrostatic pressures, and (2) the underlying physical mechanisms responsible 
for this overlooked phenomenon. 

II. COMPUTATIONAL DETAILS. 



The ah initio calculations are performed in the framework of the density functional 
theory, mostly within the local-density approxima tion (LDA) but a l so sometimes within 



1986 



the g radient generalized approximation (GGA) ([Perdew and Yue 
1996). We have employed the pseudopot e ntial method using the ABINIT 



Perdew et al. 



2002), CUSP JBel 



PWscf fBaroni et al 



aiche 



and Vanderbilt 



2000; 



King-Smith and Vanderbilt . 



Gonzc 



et al. 



mm, and 



20011 ) packages. The calcu 



norm-conserving pseudopotentials of Teter (jTeter 



ations have been done with the extended 



1993) or with the Vanderbilt-type ultra- 
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soft pseudopotentials (|Vanderbiltl . Il990r ) - always incorporating the semi-core states in the 
valence, unless specified otherwise. The wavefunction has been expanded in planewaves up to 
a kinetic energy cutoff of 120 Rydberg when using the extended norm conserving pseudopo- 
tentials and of 40 Rydberg whe n using the ultrasoft pseudo potentials. We use the Ceperley- 
Alder ex change and correlation ( Ceoerlev and Aider! . Il980r ) as parameteriz ed by Perdew anc 



Zunger (IPerdew an d Zun gerl Il981l) or the Teter-Pade parametrization ((Goedecker et al. 



1996 



Perdew and Wane , 



1992). For the ground state calculations, integrals over the simple 



cubic Brillouin zone have been rep 



point mesh (jMonkhorst and Pack 



aced by sums on a 6 x 6 x 6 Monkhorst-Pack special k- 



1976). We determined: (i) the total energy for different 



volumes of the unit cell allowing for full relaxation of the cell shape and the positions of 
the atoms inside the unit cell, in order to determine the evolution of some phases as a func- 
tion of hydrostatic pressure, (ii) the phonon spectrum of specific structures by using both 
finite-difference (frozen-phonon) and linear-response approaches (|Gonze et all Il992l ). (iii) 
the Born effective charges, the dielectric constants, and (iv) the electronic band-structure 
and density of states. It should be noted that we checked that our results are qualita- 
tively independent on technical details such as the used exchange-correlation functional, 
pseudopotential types, kinetic energy cutoff or softwares. 



III. RESULTS AND DISCUSSION 



A. PbTiO-3 under pressure 



Lead titanate PbTiC>3 is one of the most important and best studied ferroelectric ma- 
terials adopting the simple-cubic perovskite structure at high temperature. It is known 
that PbTiC>3 undergoes a single ferroelectric transition from the cubic to a tetragonal phase 
below 766 K. In fact, it appears to be a textbook exa mple of a displacive ferroelectric 



transition (Burns and Scott 



1973 



Lines and Glass . 



19771 ). 



Hydrostatic pressure P r educes, and even annihilates for high e nough value, ferroelec- 



tricity (FE) in perovskites ()Sani et al 



2004; 



Saniurio et al. 



1983), according to the con- 



ventiona l theory. Such annihilation seem s consistent with experimental X-ray diffraction 



studies (Sani et al 



2002: 



Zha et al. 



19921 ) and first-principles calculations ()Wu and Cohen , 



2005T ) that reported a ferroelectric to paraelectric phase transition in lead titanate under a 
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"m oderate" (i.e. around 10-20 GPa) pressure. On the other hand, one of our recent stud- 



ies (Kornev et al. 



20051 ) - that combines experimental techniques and ab-initio simulations 



- revealed that perovskites enhance their ferroelectricity as pressure increases above a critical 
value (which is around 30 GPa in PbTi03). 

More specifically, we performed K calculations on PbTiOs within LDA and for a wide 
pressu re range. T he calculated pressure-volume data were fitted to a Birch equation of 
state ( Birch . 198fiL Here, we present the results for 7 phases: the paraelectric cubic Pm3m 



state; the ferroelectric, rhombohedral R3m, and tetragonal P4mm phases; the antiferrodis- 
tortive tetragonal 14 /mem and rhombohedral R3c phases; and the tetragonal 1 4cm and 
rhombohedral R3c phases that can exhibit both polar and antiferrodistortive degrees of 
freedom. 

Figure HJs shows the predicted pressure dependency of the difference AH = H - H cu bic 
between the H C uUc enthalpy of the paraelectric cubic state and the H enthalpies of the 
other 6 phases. For the lowest (P) pressures, only ferroelectric distortions exist in 14cm 
and R3c, which are thus identical to PAmm and R3m, respectively. The P4mm phase 
has the lowest AH, as consistent with the fact that it is the well-known ground state of 



PbTi Os . Figs.[Tj6 ; c) show that, as expected in perovskites ((Samara et al 



1975 



Sani et all 



20041 1. significant antiferrodistortive deformations occur in the 1 4cm and R3c phases when 



increasing P above ~3 GPa. The reason leading to the TiOq octahedra rotation is clear: 
due to the TiO e octahedra rotation, the Ti-0 bonds perpendicular to the rotational axis can 
relax towards their equilibrium lengths, which minimizes the bond compression. In other 
words, the rotation relaxes the compression of the Ti-0 bond perpendicular to the rotational 
axis. On the other hand, Fig. ^ also indicates an unexpected result, namely that increasing 
pressure above a critical value enhances ferroelectricity] Indeed, an increase in P from to 
~20 GPa leads to a AH becoming less negative in th e polar phases - as consistent with the 



commonly-believed p ressure-induced reduction of FE ((Samara et al 



Sa,niurio et al 



1975 



Sani et al 



2004; 



19831 ). but the AH -vs.-P curve of these phases adopts an opposite behavior 



above J J C ~30 GPa - revealing that PbTiOs wants to become more and more ferroelectric 
after this critic al pressure. We foun d that rhombohedral states are the most stable ones at 
high pressure ((Kornev et all 120051 ) and at OK. The non-monotonic behavior of FE versus 
pressure is emphasized by the pressure dependency of the calculated c/a axial ratio in P4mm 
(see Fig. |2J): this ratio first strongly decreases with P until reaching a cubic-like value ~1 in 
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the vicinity of P c , and then it gradually increases as P further increases. Our predictions 
are thus in sharp contrast to the common expectation that the application of hydrostatic 
pressure favors a cubic paraelectric crystal structure. 



B. Possible driving mechanism (s) for the high-pressure ferroelectric instability in 
PbTi0 3 . 

In order to identify the driving mechanism(s) of the ferroelectric instability under high 
pressure, we decided to investigate in more details (1) the cubic paraelectric Pm3m phase 
and (2) the tetragonal ferroelectric PAmm phase. Fig. ^ shows that the magnitude of the 
polarization above the turning point is smaller than that below ~ lOGPa. Based on this 
fact, one might assume that the polarization is no longer the primary order parameter at 
high pressure or, equivalently, that the ferroelectric instability is driven by the coupling with 
other degrees of freedom, such as (i) the antiferrodistortive and (ii) elastic variables. 

Possibility (i) is a primary suspect because antiferrodis tortive degrees o f freedom are 
known to increase with pressure in ferroelectric perovskites ()Sani et all 120041 ) . In this case, 
the ferroelectric instability could appear as a result of a specific interaction between the an- 
tiferrodistortive and fer roelectric distortions. Howeve r, these two degrees of freedom usually 



compete in perovskites ((Vanderbilt and Zhonel 1199 8) . Furthermore, as shown in Figure Q 



the ferroelectric instability exists even without taking into account the antiferrodistortive 
degrees of freedom. For instance, after the turning point around ~ 30GPa, the P4mm phase 
(that does not contain any rotation of the oxygen octahedra) is more stable than the cubic 
phase. These facts thus rule out the antiferrodistortive degrees of freedom to be the driven 
force for the ferroelectric instability. 

Possibility (ii) is also a candidate for explaining high-pressure FE because perovskite 
ferroelectrics are well-known to exhibit a coupling between polarization and strain. An 
elastic instability may thus lead to the occurrence of a nonzero polarization. To check such 
hypothesis, we performed ab initio calculations of the elastic tensor in the cubic phase using 
the ABINIT package and found no elastic instability in the entire pressure range. Fig. El 
shows the pressure dependence of the elastic stiffness coefficients cyy , C12 , and C44. The en 
and C12 coefficients increase linearly with increasing pressure. On the other hand, the C44 
coefficient starts to decrease gradually above the critical pressure but remains positive at 
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all investigated pressures (that is, up to ~ 150 GPa) . Furthermore, these elastic stiffness 
coefficients satisfy the generalize d elastic stabi l ity cri teria for cubic crystals under hydrostatic 
pressure [see, for example, Ref. (|Karki et all Il997l )]: Cn + 2ci 2 > 0, Cn — Ci 2 > 0, c 44 > 0. 
Moreover, we found that the P4mm phase without lattice relaxation is more stable at high 
pressure than the cubic phase. Therefore, we may now assume that the polarization is 
the primary order parameter even after the turning point. Elucidating the mechanisms of 
the high-pressure ferroelectric instability should thus involve the consideration of the soft 
transverse optical mode alone. (Note that we can also conclude that the electromechanical 



coupling that beco me enormous in t 
phases of PbTi0 3 (jCohen and Wu . 



le tra nsition regions from the tetragonal to high pressure 



2006) is not responsible for the ferroelectric instability). 



C. Lattice Dynamical Properties of PbTiOs and of Related Perovskites. 

It is well-known that the Born effective charges are of primary interest for the soft TO 
mode behavior when a competition between the Coulomb and non-classical short-range 
forces is important. In the following, we will denote On the oxygen atoms located between 
two B-atoms along the z-direction and O^ the oxygen atoms located between two B-atoms 
in the perpendicular directions - when slightly displacing the atoms along the z-axis when 
computing the Born effective charges. Our calculations of the Born effective charges as a 
function of pressure in the cubic phase (see Fig. EJ) show that the anomalously large values of 
Ti and Oy persist throughout the pressure range of interest, although these effective charges 
decrease with increasing pressure For instance, their values are reduced by ~ 1.4% at ~ 25 
GPa with respect to their values at OGPa]. On the other hand, the effective charges of Pb 
and O^ ions increase with pressure, and deviate from their values at z ero pressure by ~ 2% 



at ~ 25 GPa. Interestingly, it was demonstrated by Ghosez et al. f)Ghosez et all |1996) 
that a change corresponding to a reduction of the order of 1% of the Ti effective charges 
is enough to suppress the ferroelectric instability in BaTiOs. From this point of view, the 
tendency of the PbTiOs system to become cubic as the pressure first increases from zero is 
quite understandable. However, since after the turning point the Ti and On effective charges 
continue to decrease with pressure, one might expect that the system would continue to 
stay cubic which is definitely not the case. Therefore, the pressure behavior of the effective 
charges or, equivalently, the strengths of the Coulomb interaction alone is not sufficient to 
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understand high-pressure FE. 

We now turn to Fig. |5]a that displays the LDA-simulated pressure behavior of the zone- 
center transverse optical (TO) mode frequency in the cubic Pm3m phase of PbTiOs. The 
square of such frequency is always negative, which indicates that PbTiOs has a ferroelectric 
instability and thus can not be cubic for any pressure at OK. More importantly, the mag- 
nitude of this square first decreases with P and then increases above the critical pressure 
(P c ) - confirming that PbTiOs has to become more and more ferroelectric above P c . To 
understand further the high-pressure ferroelectric instability it is advantageous to divide 
the dynamical matrix into (1) the Coulomb part (following the formalism for the rigid ion 
model ) that can easily be found once the Born effectiv e charges and dielectric tensor of 
the system have been calculated (Giannozzi et all Il99lh . and (2) the (remaining) second 
part - which is deduced by subtracting the first part from the total dynamical matrix. This 
second part thus differentiates the real material from its rigid ions picture, and gathers the 
(non-classical) short-range interactions. Figure Ef) shows the calculated contribution of these 
two parts to the zone-center soft TO phonon frequency in the cu bic phase of PbTiOs and 
BaTiOs, as a function of P. Below P c and as commonly expected ([Samara et all I1975T ). the 
Coulomb interactions favor FE while short-range interactions tend to annihilate it - since 
they have a negative and positive frequency square, respectively - and the decrease of FE 
when increasing P is driven by the cost of short-range interactions becoming predominant 
over the gain associated with long-range interactions. Figure Ef) also reveals that FE above 
P c dramatically differs from the low-pressure "normal" ferroelectricity because of the reverse 
sign of Coulomb and short-range parts: it is now short-range effects (which are electronic 
in nature) that favor FE at high pressure while Coulomb interactions (that involve ions 
"dressed by electrons") prefer paraelectricity! 

Figure^ indicates that the results that we present are qualitatively general for perovskite 
ferroelectrics and not just an artifact of studying lead-contained perovskites, since they also 
apply to BaTiOs. This fact is further emphasized in Figures |B1 that report our OK ab-initio 
predictions for the eigenvalue A and eigenvectors (£4, £b, £p,, £ , £o» ) of the second- 



derivative matrix defined in Ref. ((King-Smith and Vanderbilt 



1994) and associated with the 



soft-mode in cubic PbTiOs, BaTiOs an d BaZrOs, as a function of pressure. (It is important 
to realize that the last two materials are lead-free). A negative A corresponds to a ferroelectric 
instability while the magnitude of the eigenvectors quantifies the contribution of the different 
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ions to FE, if any. All the data displayed in Figs. El are obtained by using the local density 
approximation (LDA) except in parts (a) and (d) for which we also report results from 
the use of the generalized gradiant approximation (GGA). Figs. El undoubtedly show that 
all studied materials should be ferroelectric at high enough pressure, thus confirming that 
one has to rule out a possible suggestion that pressure-induced FE only occurs in Pb-based 
materials. One can further notice that £4 has a very small magnitude at high pressure 
in PbTiO"3, BaTiOs and BaZrC>3. This fact clearly implies that the A atoms contribute 
very little to high-pressure ferroelectricity. Our Figs. El also nicely summarize the different 
behaviors one can get as pressure increases in transition-metal perovskites, as a result of our 
discovered non-monotonic pressure behavior: the existence of FE at any pressure (PbTi03, 
see LDA calculation of Fig. EJa)), the disappearance and then reentrance of FE under 
pressure (BaTiOs, see Fig Efb)), and the occurrence of FE at high pressure in a nominally 
paraelectric compound (BaZrOa, see Fig. Etc)). Moreover, Figures Efd-f) clearly shows 
that the eigenvectors of B and O atoms have a large (and mostly material-independent!) 
magnitude at high pressure, which points to an "universal" microscopic explanation for 
high-pressure FE. Finally, such pressure-induced phenomenon is also found when performing 
ab-initio calculations in KNb03, NaNb03, SrTi03, CaTi03, as shown in Fig.0 Pb(Zr,Ti)03 
and some Bi-based perovskites (not shown here). Not also that the fact that the A atom 
contributes little to high-pressure FE is clearly confirmed in Fig. [SJ-that shows that A-free 
WO3 material also enhances its FE at high enough pressure. 

D. Electronic Structures under Pressure 

We now resume our "detective" work and look for the microscopic origin of the observed 
high-pressure FE. As it was shown above, the instability is related to the non-classical short- 
range interactions. As a result, we will now mostly concentrate on electronic properties. 

1. Influence of semicore states. 

We first investigate the role played by the semicore states via Fig. that shows the pre- 
dicted pressure behavior of phonons in PbTiOs in the cubic phase when neglecting to incor- 
porate the semicore states of Pb and Ti as valence electrons. Technically, the ab-initio sim- 
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ulations without the semicore states have been performed with the Troullier-Martins norm- 
conserving pseudopotentials within ABINIT. Fig. ^indicates that the semicore states Pb 5d w 
and Ti 3s 2 3p 6 (1) must be included explicitly to produce the well-know n ferroelectric instabil- 



ity at zero pressure - as consistent with p revious first-principles works ((Cohen and Krakauer . 



1992: 



King-Smith and Vanderbiltl . ll994l ) but (2) are not the reason for the existence of a turn- 



ing point and of high-pressure FE (Note, however, that the magnitude of the instability is 
of course not the same when neglecting or incorporating these semicore states in the valence 
electrons). Note that the fact that Pb 5<i 10 orbitals are not responsible for the existence of 
the ferroelectric instability at high pressure could also have been guessed thanks to Fig. |H1 
that shows that tungsten trioxyde WO3 (which is a defect-perovskite structure having no 
A- atom) also exhibits a highy-pressure FE. In fact, this latter figure, as well as Figs. (|6I7|) . 
indicates that the electrons of interest for high-pressure FE should belong to the O and/or 
B ions. 



2. Influence of the d-states of B atoms. 

Figure El reports the LDA-predicted pressure behavior of the soft-mode eigenvalues of 
the second-derivative matrix in cubic bismuth aluminate BiA103 and magnesium silicate 
MgSi0 3 . Note that both materials, unlike the ones considered above, do not possess any 
B-atom having occupied d electrons. Interestingly, no turning point is observed in MgSiOs - 
under pressure up to 300 GPa. As a matter of fact, we found that the TO soft mode at the 
T point linearly increases with pressure in this material and become positives at around 3 



GPa - which is consistent with previous studies of MgSiO"3 ([Parlinski and Kawazoe . 2000). 
The same qualitative behavior is found in BiA10"3 under high pressure up to 900 GPa. We 
thus arrive at the conclusion that the d states of the B atoms are of primary importance for 
understanding the physics behind the high-pressure ferroelectric instability. 



3. Electronic band structure and density of states under pressure 

To further elucidate the nature of the high-pressure ferroelectric instability and the role 
played by the electronic structure, we have calculated the total and atomic site projected 
densities of states (PDOS) of ideal cubic perovskites. Figs. [TT] display some partial electronic 
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density of states for cubic PbTiOs under different pressures. Note that the contribution 
of Pb orbitals is not shown because, as discussed above, Pb was found to be relatively 
ferroelectric-inactive at high pressure. 

Figures UTT q-c) show the low-lying states ranging between -22 and -15 eV (that solely 
originate from the O 2s orbitals in an ionic picture), as well as the topmost group of valence 



bands (that exhibits the we l 



199 



2: 



Cohen and Krakauer 



-know n hybridization between O 2p and Ti 3d orbitals ((Cohen , 



1990)). As it is known, the topmost states correspond mainly to 



a (bonding) and 7r (non-bonding) orbitals. Moreover, the states with large Ti 3c? contribution 
are mainly concentrated in the lower (bonding) part of the topmost region. This O 2p and Ti 
3d hybridization persists throughout the pressure range of interest. The integrated PDOS of 
this group have a smooth and small pressure-dependency thus ruling out that high-pressure 
FE is caused by this p - d hybridization. At the same time the Ti 3d state contribution to 
the low-lying O 2s-\ike bands has increased considerably. Thus, as the pressure increases, the 
admixture of Ti 3d states in the occupied bands seems to evolve differently for the occupied 
low-lying and topmost bands. More precisely, the evolution of O 2s - Ti 3d and O 2p - Ti 3d 
mixing is nicely seen in the trends displayed by the integrated (over all the occupied valence 
bands) PDOS (to be denoted by IPDOS), see Fig. H2». The IPDOS of O 2s and Ti 3d 
vary in a correlated way with pressure, the former decreasing and the latter increasing. 
On the other hand, the IPDOS of O 2p remains more or less constant, apart from the 
highest pressures. These tendencies are even more clearly seen in Fig. I12f ft), where IPDOS 
are normalized with respect to their corresponding values at zero pressure. Furthermore, 
Fig. lT2T c) reports the contributions to the IPDOS coming from the low-lying states (located 
between -22 and -15 eV) and from the topmost group of the valence bands. The most striking 
effect is that Ti 3d orbitals (and more precisely their e g parts) dramatically increase (in a 
nonlinear way) their mixing with the O 2s orbitals above P c . The admixture of Ti 3c? states 
in the occupied levels of the electronic structure releases O 2s - O 2s antibonding states 
above the Fermi level, thus optimizing the band energy. A similar behavior of integrated 
PDOS was found in BaTiOs under pressure, as reported in Fig. [THl 

A rapid relative increase of the O-s and B-d mixing can easily be understood as follows. 
The lobes of the Ti 3c? — t 2g orbitals point between the oxygen ligands (whereas the Ti 
3c? — e g orbitals point directly towards the ligands), forming the a bonding. Hence, the 
overlap with O 2s orbitals will be significant at high enough pressure for the Ti 3c? — e„ 
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states, which results in local repulsion between overlapping charge densities. From one 
hand, this repulsive interaction will cost some energy. From the other hand, some of the 
repulsion is compensated by hybridization in the resulting bonding states - which in turn 
leads to considerable amount of O 2s state in the conduction band. Taking into account 
the different interatomic dependence of these two competing interactions, one might expect 
that the need to avoid the large local repulsion will win at some pressure, enabling atomic 
distortions and thus giving rise to ferroelectricity. 

We would like to point out that the important role played by the small amount of 2s 
admixture with the e g orbitals of transition-metal cations surrounded b y anions (like chlorine , 



bromine, fluorine, or oxygen) was a. 



Moreover, recent DFT calculations ((Hussermann 



so previously recog nized in Ref. ([Moreno et all \1997). 



20031 ) revealed the key role of s — d mixing 



on the stability of high-pressure structures in P, As, Sb, and Bi. 



E. Topological analysis of the electron density in PbTiC>3. 

To get more insight into the chemical bondings of PbTiC>3 under pressure, we consider the 
topological properties of the electron density p(r) . The total density was restored by adding 



atomic core electron density to the valence density obtained in pseudopotential calcu 



The analysis was performed in terms of the atom-in-molecules Bader theory (jBaderl . [1990) 



at ions. 



which is based upon the so-called critical points, r cp , of the electron density. At these 
points, the gradient of the electronic density is null and the Hessian matrix d 2 p(r) / 'dxidxj 
is characterized by its three eigenvalues, Interatomic interactions can be adequately 
described by the topological properties of the electron density at the so-called (3,-1) critical 
points (or bond critical points) for which the three eigenvalues of the Hessian matrix are 
non-zero and the algebraic sum of their signs is —1. Bond critical points occur as saddle 
points along bond paths between pairs of bonded atoms and can be used to define atomic 
coordination. The characteristics of the bond critical points in PbTiOs at GPa are listed 
in Table 1. 



As expected for cubic perovskite [see, for instance, Ref. (jZhurova and Tsirelsonl . 120021) ]. 
we found six (3,-1) critical points on the Ti-0 lines and twelve (3,-1) critical points on the 
Pb-0 lines at any pressure. Interestingly, all the (3, -1) critical points in PbTiOs are charac 



terized by the positive Laplacian of the electron density which, according to Bader ((Bader . 
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TABLE I Characteristics of some critical points (CP) in PbTi03 at a = 7.35 a.u. (zero pressure); 
V 2 p(r cp ) = Ai + A2 + A3, where Ai and A2 are perpendicular to the bond line curvatures of the 
electron density and A3 is parallel. Values are given in atomic units. 



Top 


p(r cp ) 


V 2 p(r- cp ) 


Ai 


A2 A3 


Type of CP 


(0.234, 1/2, l/2) a 


0.132 


0.590 


-0.219 


-0.219 1.028 


(3,-1) 


(0, 0.267, 0.267) 6 


0.028 


0.078 


-0.020 


-0.022 0.120 


(3,-1) 



a Ti-0 line. 
^Pb-O line. 



199(1 ). is an indication of the closed-shell interaction between the corresponding atoms. 



Analyzing the topologica l characteristics available in the literature, it was found 



( Zhurova and Tsirelson 



2002) that the bonds, typically identified as ionic, are characterized 
by the electron density value at the (3,-1) critical point in the range of 0.07-0.25 e A -3 and 
a specific range of the ratio of the perpendicular and parallel curvatures of the electron 
density: 0.12 < |Ai|/A 3 < 0.17. According to Table I, the Pb-0 bond (at GPa) fits well 
the first condition since p{r cp ) m 0.1908 e A -3 and closely fits the second one if we take 
|(Ai + A 2 )/2|/A 3 « 0.176 rather than |Ai|/A 3 . Thus, the Pb-0 interaction (at zero pressure) 
can be considered as mostly ionic with small covalent admixture. On the other hand, the 
value of p{r cp ) ps 0.9 e A -3 and the ratio | Ai | / A3 = 0.213 at the (3,-1) critical points on the 
Ti-0 line imply that the Ti-0 bond is not ionic in nature (at zero pressure). 

In order to characterize the chemical bonds in more details (and for different pres- 



sures), we emp 



(Es pinosa et al. 



oy the recently proposed cla ssificat i on sc heme of interatomic interactions 
2002 ). According to Bader ( Baderl 1990 ) the Laplacian function V 2 p(r c 



cpj 



is related to the potential energy density v(r cp ) (of negative value) - which represents the 
capacity of the system to concentrate electrons at the critical point - and the kinetic energy 
density g{r cp ) (of po sitive value) - which g ives the tendency of the system to dilute electrons 



at the same point (jEspinosa et al 



1998) - through the local form of the virial theorem 



l/4V 2 p(r cp ) = 2g(r cp ) +v(r cp ). Note that for pressure different from zero, the virial theorem 
should include an additional term proportional to this P pressure: in the cubic case, one 
should add the +3P term to the lefthand side of this latter expression. 

The kinetic energy can be calculated using the following expression for closed-shell in- 
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teractions at the bond critical points (|A bramov . 



g(r cp ) = 3/10(37T 2 ) 2 /3 p ( rcp )5/3 + 



l/6V 2 p(r cp ). The potential energy density v{r cp ) is then calculated by subtracting g(r cp ) 
from l/4V 2 p(r cp ). We also computed the total electronic energy density h(r) at the 
r cp points as h{r cp ) = g{r cp ) + v{r cp ). Using the adimensional \v(r cp )\/ g(r cp ) ratio and 
the bon d degree BD = h(r C P ) / p(r rp ), the atomic interactions can be divided into three 
classes (E spinosa et oil Il998l ): (i) pure closed shell for which \v{r cp )\/ g{r cp ) < 1 (and for 
which the V 2 p(r cp ) > and h{r cp ) > inequalities stand), (ii) pure shared shell for which 
1 < \v{r cp )\/ g{r cp ) < 2 (and for which V 2 p(r cp ) < and h{r cp ) < 0) and (hi) transit closed 
shell for which \v{r cp )\/ g{r cp ) > 2 (with V 2 p(r cp ) > and h{r cp ) < 0). 

Inside the ionic region (i), the bond degree parameter is an index of ionicity: the stronger 
the ionicity the greater the BD magnitude. In regions (ii) and (iii), the BD parameter 
measures the covalency degree: the stronger the covalency the greater the BD magnitude. 
For instance, the ionic Li-F bond is characterized by \v(r cp )\/ g(r cp ) pa 0.76 and BD pa +0.41 
while t he covalent C-0 bond in CO(NH2)2 is such as \v(r cp )\/ g(r cp ) ~ 2.64 and BD pa 



—0.29 (Tsirelson, 



2002). Fig. fT4l reports the bond degree parameter versus \v(r cp )\/ g(r 



cp I 



ratio (as determined from our hrst-principles calculations) for the Ti-0 and Pb-0 bonds. 
This figure confirms that the Pb-0 bond lies close to the pure closed shell region (i.e., it is 
mostly ionic) at OGPa, and reveals that such bond becomes more ionic at higher pressure. 
On the other hand, the Ti-0 bond lies close to the middle of a transit region (ii), implying 
that the Ti-0 bond is neither purely covalent nor purely ionic. Fig. El also shows that 
the covalency degree of the Ti-0 bond increases with pressure, which is consistent with 
our band structure calculations discussed above in subsection D.3. Interestingly, since this 
covalency increases with pressure, there might be a possibility for the phase transition to 
occur due to the drastic collapse of the oxygen ionic size (which is related to charg e transfer 
to oxygen). For example, the Shannon ionic rad ius of 2 ~ is 1.23 A f Shannon! . 1 1 9 76 ) . whereas 
its covalent radius is 0.73 A ( Sanderson! . Il962 ) . If the size of oxygens were small compared 
to the allowed space then the oxygen atoms would be weakl y bounded in the lattice and 
a multi-well pot ential could be expect ed [similar to PbTe:Ge fjMurase , 



and PbTei-x^a, rtAbdullin et al 



1980; 



Weiser 



1981) 



1984)]. The mechanism for the off-center instability of the 



O ion could be viewed in this case as an imbalance between the decreased (due to small 
O radius) repulsive forces and the polarization forces that tend to displace the ion from 
its position. Furthermore, the on-site force constant of the O ion would become negative 
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and the O ion would be unstable to the displacement from its regular high symmetry site 
with local dist ortions. This would le ad to the instability of O similar to one observed for 



Li in LiN bO.q (Parlins 



volumes ( Bader , 



1990; 



ri et al 



Gonze et al 



2000) . Ind eed, although we found that the so-called Bader 



2002) (which are a measure of the ionic volumes) are 



smooth functions over the pressure interval of interest, the Bader O volume, unlike those 
of Pb and Ti, decreases faster under pressure than the volume of the unit cell, as shown 
in Fig. However, our ab initio lattice vibration calculations show that the on-site force 
constant of O is positive at any pressures and, hence, the local atomic potential possesses a 
single minimum. We thus can rule out the drastic collapse of the oxygen ionic size as the 
mechanism for the high-pressure ferroelectricity. 



F. Simple Tight-Binding Model 

To confirm the role of the interaction between the 2s orbitals of oxygen atoms and the d 



orbitals of the B atoms at high pressure, we use the bond-orbital model of Harrison ({Harrison . 



19801 1 - a simplified tight-binding model, where the Hamiltonian is limited to the on-site and 
nearest-neighbor terms. This approximation allows an extraordinary simplification of the 
problem and the total energy can be rigorously calculated bond by bond. This is a very crude 
approximation but, surprisingly, it works very well in a wide variety of systems. In particular 



it exp lains reasonably well the dynamical charges in ferroelectric perovskites (jGhosez et al 
19981 1 . In this section, we shall make use of the method to clarify the role of the O 2s - B d 
and O 2p - B d couplings in the high pressure ferroelectric instability in perovskites. We will 
see that the instability at high pressure is governed by a competition between short-range 
covalent and overlap interactions. 

Intuitively, the gradually decreasing ferroelectricity with rising hydrostatic pressure from 
zero to small values (i.e., before the turning point) may be attributed to the fact that 
the conventional mechanism responsible for triggering the ferroelectric instability at zero 
pressure and related to the B d - O 2p hybridization is not working well under pressure 
because of an increase of covalency. It is known that in the case of significantly covalent 
systems a paraelectric phase is of lower energy than the ferroelectr ic one due to the increasin g 



of the rigidity and the stability of the metal-oxygen network ([Villesuzanne et al 



19981 1 



From this point of view, it seems likely that - besides the pd covalency that appears to 
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prevent the ferroelectric instability at high pressure - there is an interaction that reveals itself 
only at high enough pressure (i.e., when the distance between atoms is sufficiently small). 
This gives us a clue that the interaction should be small at large interatomic distances (e.g., 
at zero pressure) and should increase more rapidly than covalency effects at small interatomic 
distances, bringing the ferroelectric instability into the system. 

For a qualitative understanding of phase transformations under pressure, we recall that 
two energy contributions to the cohesive energy, namely the covalent and overlap energies, 
depend most strongly on the interatomic distance. According to the Harrison model, the 
(attractive) covalent interaction depends on the interatomic distance, d, as ~ d~ 7 ^ 2 while 
the (repulsive) overlap interaction - which is due to nonorthogonality of wave functions on 
adjacent sites - behaves as ~ d~ 7 and thus depends more strongly on the nearest-neighbor 
distance. The interplay between the short-range covalent and the overlap interactions can 
be nicely seen from the following arguments. The covalent energy contribution mainly 
represents the short-range electron-ion potential energy while the overlap energy can be 
predominantl y ascribed to th e increase of the kinetic energy of the valence electrons upon 
compression ( Harrison . 198(tL Since the resulting forces due to these two interactions (cen- 



tral repulsive forces and covalent forces) scale differently under pressure, one might expect 
that at some pressure the noncentral nature of the bonding will become negligible. Indeed, 
Fig. EI] shows that the high-pressure Cauchy relation between the elastic constants for hy- 
drostatic conditions C12 = C44 + 2P, where Cu and C44 are the elastic constants calculated 
from first-principles, is satisfied for cubic PbTiC>3 in the vicinity of the the turning point 
(~ 30GPa), which indicates an insignificant role of noncentral and many-body forces in this 
pressure region. This fact implies that, for instance, empirical two-body potentials could be 
used for calculations of elastic properties for pressures around the turning point. Away from 
the critical pressure the deviation from the Cauchy relation becomes significant. The latter 
implies that the contribution of either noncentral or many-body forces becomes more and 
more important at higher pressures, and it cannot be treated anymore as a small correction 
to the two-body potentials. Thus, we clearly see manifestation of the competing forces under 
pressure. 

Note that a major role of competi tion between covalent and overlap interactions has 



been recognized for sp-bonded solids ([Maiewski and Yogi 



19851 ). Remarkably, a drastic 



softening of the transverse optical phonons across the pressure-induced phase transition 
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was predicted, and the physical origin of this softening was shown to be closely related to 



ferroelectricity ([Maiewski and Vogl . 



1987). 



Let us now proceed further by approximating, in Ti-based perovskites, the total effect 
of the coupling (i) between the oxygen 2s states and the d-orbitals of the Ti atom, and (ii) 
between the oxygen 2p states and the <i-orbitals of the Ti atom, via the following expressions: 



E sd bond = 2V 3 (sd) - 2{V 2 {sd) 2 + Vsisd) 2 ] 1 / 2 + E sd rep 
E pd bond = 6V 3 {pd) - 6[V 2 {pd) 2 + V 3 (pd) 2 f 2 + E pd rep 



(1) 



where the (polar energies) V 3 are the average of the cation and anion energy, and are 
given by V 3 (sd) = (e d - e s )/2 and V 3 (pd) = (e d - e p )/2 - with e s = -29.14 eV, e p = -14.13 
eV, and e d = —11.04 eV (which are the energies of the 2s and 2y states o f the oxygen 
atom and of the 3d states of the titanium atom, respectively ((Harrisonl [l980j). The second 
terms f...] 1 / 2 (which are the covalent energies) are connected with the bond formation energy 
and are responsible for l owering the energy of the bonding state. Based upon a second- 
moment approximation ((Wills and Harrison! . 119831 ) . these covalent energies are given by 
V 2 (sd) 2 = EjVsd^rj) 2 and V 2 (pd) 2 = T, j {V pda (r j ) 2 + 2V pd7V (r j ) 2 )/3, where the last sum is 
over the 6 nearest neighbors to Ti and where the so-ca lled Harrison's two-center interactions, 



Vidm, can be taken to be of the form ( Harrison 



1980): 



V, 



Idm 



Vidm 7 in 

V rV 



(2) 



where rj is the interatomic distance between Ti and the j th nearest-neighbor O, fi is the 
free-electron mass, r d = 1.029 A for titaniu m, m = cr.7r, and r/'s take the following values: 



Vsda = -3.16, r] sdu = -2.95, and r] sda = 1.36 ((Harrisonl. Il980h. 

The third terms of Eq. (1), E sd rep and E pd rep , are related to the overlap interactions. 
A repulsion comes from the shift in the energy of each electron due to the nonorthogo- 
nality of the s, d and p, d s tates on adjacent sites, with the average shifts being given by 



(Harrison and Straub 



19871 ): 



Se sd = -^2V sda (rj)S sda (rj) 
j 

5e pd = -1/3^ \Y pdtT (rj)S pdu {rj) + 2V pd7T (r j )S pd7r (r j )] 



(3) 
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The overlap repulsion for a given bond, is equal to this shift times the number of elec- 
trons i n these sta t es. I n the approximations adopted here we use the extended Hiickel 



theory f Ho ffmann!. 



19631) to derive the overlap eleme nts, as proposed by van Schilfgaardo 



and Harrison ((van Schilfgaarde and Harriso 



1986). It is plausible to take the overlap 



elements to be proportional to the corresponding off-diagonal hopping elements and in- 
versely proportional to the on-site atomic orbital energies: S a d m (r) = ■j^—V a dm(r)/(e a +eA 

-* *- ad 

where a = s,p and m = a, n. For the Wolfsberg-Helmholz parameter K nd, the weightec 



formu la K ad = k + (e a — ed) 2 /(e Q + e^) 2 + (1 — K)(e a — ed) 4 /(e a + e^) 4 (jAmmeter et al. 



19781 ) could be used. In gene ral, the Wolfsberg-Helmhol z parameter K a d can be given by 



the distance-de pendent form (ICalzafe rri and Rvtz 



following Refs. ( Baranowski . 



198 



Harrison! . [1983; 



199fiL 



For the sak e of si mplicity and 



Maiewski and Vogl 



1987), we assume 



that K takes a different value for each row of the Periodic Table, i. e., the K of Ca rbon 
(Kq = 1.63) for Oxygen and the K of Ge (Kq £ = 1.07) for Titanium ((Harrison! . Il999j ) and 
take the average as K = (KcKce) 1 ^ 2 - 

Following this scheme, the covalent or polar bonding character of a given bond (e.g., sd 
or pd) can now be addressed. For each bond, the bond polarity, a p , of the compound can be 
obtained from a 2 , = V$ /(V 2 + V 2 ) and a complementary quantity, the covalency coefficient, 



ot c , is defined through the relation a c = 1 — a^. The tight-binding model demonstrates the 
strong reduction of the polarity of the sd and pd bonds under pressure, as shown in Fig. El 
On the other hand, the model predict a considerable increase of the sd O — Ti chemical 
bond covalency under pressure (i.e., when the interatomic distance decreases) - which is 
consistent with the first-principles population analysis reported in Section III.D.3. Indeed, 
as it can be clearly seen from Fig. El the sd bond covalency depends more strongly on 
the nearest-neighbor distance than the pd bond covalency, and therefore the first one may 
predominantly determine the behavior of perovskites in the vicinity of the turning point. 
In the Harrison bond-orbital for malism, the frequency of the q = optical phonon can 



be easily obtained (H arrison! . Il980f ). Consider a perovskite structure viewed along a (100) 
direction with a relative displacement Uj = un — uoj of the Ti and O atoms from their 
equilibrium positions in a (100) direction when the atoms oscillate according to a T-point 
soft mode. Accordingly, the change in length of the indicated bond Ti — O is 5dj, leading to a 



change in the bonding energy of SE sd bond and SE } 



pd bond' 



These changes in the bonding energy 



may be viewed as the frozen-phonon energy obtained as the difference of the respective energy 
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relative to its equilibrium value, 5E ad bond = E ad bond(uj) — E ad bond(0). For the evaluation 
of the sd and pd contributions to the frozen-phonon energy, usual harmonic fit function can 
be used, i.e. 2 5E ad bond = uS^ d J2j Mju 2 , where Mj are the atomic masses, and u> ad are the 
contributions from sd and pd interactions to the total frequency. [As an alternative to the 
above method, an analytic expression for these changes can be obtained by expanding Eq. Q 
in a Taylor series with respect to the displacements and then truncating the series to the 
second order in the displacements.] 

Let us consider the frequency of the q = optical mode in which the Ti and O atoms move 
out of phase with the relative displacements which reproduce a high-pressure eigenvector of 
the soft mode obtained from first-principles (see, e.g., Fig.6). Since above the turning point 
the relative displacement of each of the 4 transverse oxygens u± is greater than that of each of 
the 2 longitudinal oxygens u\\, 2u± > u\\, we have chosen the following displacement pattern 
along the z axis for the 6 oxygen ions: [0Au,0Au,u,u,u,u] with u = 0.002. Denoting d as 
half of the cubic lattice constant We now see that one neighbor distance becomes d — 0.4u 
(between the Ti atom and the top O atom), one becomes d + OAu (between the Ti atom and 
the bottom O atom) and four become (d 2 + u 2 ) 1 ^ 2 (between the Ti atom and 4 transverse 
O atoms). 

The results of the frozen-phonon calculations are shown in Fig. At large interatomic 
distances both sd and pd contributions stabilize the paraelectric state, since they both 
have a positive SE. The curves of SE sd bond and 5E pd bond have maximums which lies at 
around 2 A and 1.8 A, respectively. One thus finds a maximum at around 1.85 A for the 
sum of two contributions. [Note that the experimental equilibrium interatomic distance in, 
for instance, SrTiC>3 is 1.95 A and the turning point is around 1.85 A according to our 
first-principles calculations.] These estimates are compatible with the assumption that the 
repulsion overlap energy does not become important until the turning point. The most 
striking result is that the model predicts a strong phonon softening upon compression since 
Fig. CHI shows that both 5E sd bond and SE pd bond, that first increase with pressure and 
reach a maximum at around the turning point, then decrease rapidly under compression 
and become negative at m 1.83 A and ~ 1.6 A, respectively. The interatomic distance at 
which the sum of both contributions becomes negative is m 1.64 A, to be compared with our 
ab-initio result of ~ 1.78 A for SrTiC>3. [It should be emphasized that an exact agreement 
with ab-initio calculations cannot be expected, on account of the approximations made in 
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the tight-binding model, especially the neglect of long-range interaction.] Interestingly, it 
is the sd contribution which first destabilizes the paraelectric phase under the influence of 
pressure. Therefore, the phonon frequency is reduced under pressure after the turning point 
as much as this destabilization is not compensated for by the other energy contributions. 
Moreover and interestingly, as the system is further compressed the pd contribution also 
contributes to the destabilizing of the parelectric phase. The present tigh-binding model 
thus predicts a striking softening of a T-point phonon mode under compression, and clearly 
explains the strong preference of perovskites to form ferroelectrically distorted structures 
at high pressure due to an electronically driven structural instability (i.e., electron-phonon 
coupling). In particular, this simple model reveals the predominant role played by the 
costly-in-energy overlap between the 2s orbitals of oxygen and the 3d orbitals of Ti on the 
high-pressure FE of perovskites. 



G. Electronically Driven Instability 



According to the above simple tight-binding model, it is the electrons and their inter- 
actions with ions (i.e., electron-phonon coupling) that are responsible for the high-pressure 
ferroelectric instability of perovskites. In this section, the electron-ion interaction will be 
considered in more details. More precisely, we shall be interested in the electronic contribu- 
tion into the ferroelectric soft mode within DFT and within the adiabatic approximation. 

Note that the adiabatic approximation consists in (a) neglecting the phonon contribu- 
tions to the electronic density response and (b) approximating the dynamical response by 



a static one (Geilikman 



1975 



Sham 



1974 ). In the adiabatic approximation, the electrons 



are considered to remain in the instantaneous (electronic) ground state when the lattice 
vibrates. Consequently, the effect of electrons on the adiabatic phonon frequencies is due 
to the fact that the change in the charge density induced by the displacements of nuclei 
screens the e xternal perturba tion (and the linear-response equations must be solved self- 
consistently) (|Savrasovl . 119961 ). 

To refresh readers' memory, let us remind some basic equation of density-functional theory 
(DFT). Within the framework of DFT, the total energy can be expressed in terms of the 
Kohn-Sham eigenvalues e^: 
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2 r , , ,n(r)n(r') 



Jn{r)v xc {r)dr + E N {K). (4) 



where fa is the occupancy of the state i, and n(r) is the density of the valence electrons. 
The first term is the band structure energy: 



EM = £/<<*! - 1/2V 2 + v SCF (v)\i) 

i i 

where an effective potential (also called self- consistent, SCF, potential) is: 



V SC F(r) = V ext (v) + e 2 



-dr' + v xc (r) 



(5) 



(6) 



and v xc (r) = 5E xc /5n(r) is the exchange-correlation potential and V ext (r) is the external 
(ionic) potential acting on the electrons. -Ejv(R) is the electrostatic interaction between 
different nuclei: 



1+3 



Note that, in general, the bare ion-ion contribution into phonon modes can stabilize as 



well as destabilize the paraelectric state. 



trivial consequence of Earnshaw's theorem (jEamshaw, 



he destabilizing nature of this contribution is a 



1842 : 



Fevnman and Sands 



1975|) that 



states that a static system of electric charges can not have a position of stable equilibrium. A 
a result, the crystals stability is ensured by the electron contribution to the force constants. 
Although it is possible that the ion-ion contribution to a ferroelectric soft mode destabilizes 
this mode, it is not necessarily the cas e - implying that the ferroelectric soft mode can be 
stable in a lattice of bare ionic charges (|Kvvatkovskii and Maksimovl . 119881 1 . This particular 
situation occurs in the case of the ferroelectric soft mode in perovskites under pressure, since 
we numerically f ound that the ionic term in the force constants (which arises from the ion- 



ion Ewald term (|Giannozzi et all Il991[ ) with the bare ionic (pseudo) charges (for example, 
Z Pb = 14, Z Ti = 12, Z = 6)) is positive in the pressure range of interest. 

Next, we are interested in the electronic contribution to the phono n spectrum. The total 



secon d-order change in the electronic energy E el can be written as (Jo hnson 
1969(). 



1969; 



Sham . 



E. 



(2) 
el 



p(2) p (2) 
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1/2 1 dr5n(r)5y ext (r) + J drn\r)8V^ t {v) 



(8) 



where n°(r) is the electronic density of the unperturbed crystal, 5n(r) is the change in 
the electron density due to the ionic displacements, 5V ext (r) is the change in the "bare" 
potential experienced by the electron system to first order in the ionic displacements, 
and (r) the second order change. Therefore, the expression for the electronic con- 
tribution to the adiabatic force constant matrix consists of a negative-definite contribution 
E^l (t he first term in Eq. El) and a positiv e-definite contribution E^l (the second term in 
Eq.ED (|Kvvatkovskii and Maksimovl . Il988l ). 



(2) 

Let us now consider the n egative-definite contribution E„,{ to the adiaba t ic for ce con 



stant matrix in more details (Geilikman 



1975 



Kvvatkovskii and Maksimov . 



19881) . This 



contribution describes how a first order change in the external Hamiltonian leads to a first 
order change in the density matrix, which acts back at the Hamiltonian. 

Let us first consider the change in density 5n(r) induced by a change in the external 
potential 5V ex t(r). Within DFT, one can define the non-interacting electron polarizability, 
Xo( r ? r ')> as the charge- density response to a variation of the total potential: 



5n(r) = J dr' x (r,r') SV S cf{^)- 



(9) 



The expression for xo{ r , r ') has the well known form (jPick et all 11970c IShaml Il969. 19741): 

Xo (r,r') = E^— ^ ^(r)^-(r)^(r')^(r') (10) 



i,j 6 * € j 



where the sums over i and j extend to both occupied and empty states. Clearly, only terms 
involving virtual transitions from occupied or partially occupied to empty or partially empty 
states contribute because of the difference fi — fj in the numerator. 

Recall the relationship between the bare and the self-consistent perturbing potential: 



SV SCF (r) = 5V ext (r) + J dr' K(r,r')5n(r'), 



where: 



K(r,r') 



+ 



5 2 E XC 



(11) 



(12) 



|r — r'| 5n(r)5n(r') ' 

Combining Eqs. (JH HTTj) . one can rewrite the the negative-definite contribution as follows: 
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E%1 = 1/2 J dr 6n(r)6V ext (r) 

= 1/2 J drdr'5V SC F(r)xo(r,r')5VscF(r') - 1/2 J dr' 5n(r)K(r,v')8n(r') (13) 



Substituting the results (|12|) and (|13jl in Eq. (jSJ) we obtain 



4 2) = 1/2 1 rfrrfr' ( 5V5 CF (r) Xo (r,r / )^5CF(r , ) + / dm°(r)<^ 



(2), 
ext ' 



- 1/2 /" dr' <Jn(i 



+ 



5n(r') 



(14) 



<5n(r)<5n(r') 

The first two terms in Eq. (fPfj) are due to the variation of the band energy term in the 
total energy (the first term in Eq. (J3|) and the third term in Eq. is due to the variation 
of the Hartree energy and the exchange-correlation energy. 

The total electronic contribution (see Eq. [SJ) to the interatomic force constant ten- 
sor between the atoms at and ry can be linked to the second-order change in 
the electronic energy with respect to collective ionic displacements in the usual way as 



2E. 



(2) 
el 



J2iv ^u'Ui,aUi',/3, where Ui >a is the displacement of the Zth ion in the a-direction 



from its equilibrium position R;. 

Taking into account that SVscF( r ) — 12i $VscF( r ) / ■ u h the negative-definite contri- 

(2) 

butio n of the band structure e nergy E el { to the dynamical matrix can be written in the 



form ( Varma and Weber 



19791 ) 



*S?(q) 



1 /k+qA' ~ /kA la* l'f3 

AT 2-~i ~ . . _ r . . fi'k+qA',kAfi , k+qA',kA 



N, 



AA'.k £ k+qA' - e kA 



(15) 



where N c the number of unit cells in the crystal, and g l k , x , kA = (k. r \'\5VscF( r ) / 
are the electron-phonon matrix-elements. The sum has to be carried out, in princi- 
ple, over all bands A, A' and all k points k, k' that can be connected via the phonon 
wave vector q. Interestingly, simplifying th e electron-phonon matrix-eleme nts within 
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19791) #kA,k'A' oc 



nonorthogonal tight-binding representation as (j Varma and Weber 
{de^x/ dklA^SkA^xx' — de^xi /dli'lA^SwA^w) - where S is the overlap and A the eigen- 
vector matrix - we clearly see the dependence of the electron-phonon interaction on the 
overlap effects. An important feature of this expression noted in Refs. rfVarma and Weber , 
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1977 



1979) is the fact that the electron-phonon coupling is large where the electron veloci- 



ties are large and therefore the density of states are small. Indeed, we numerically checked 
these facts in our first principles calculations. 

Inte restingly, Eq. (Toj) is in fact the electron-phonon coupling in the so-called Frohlich 
model ( Frohlich . Il954 ) - which thus proves that such coupling is naturally included in th e 



adiabatic approximation within DFT (|Geilikmad . ll975clKvvatkovskii and Maksimovl .ll988). 



(Note that the electron-phonon cou pling is responsible for the fe rroelectric instability in the 



vibronic theory of ferroelectricity ([Kristoffel and Konsin . 
pseudo Jahn- Teller effect due to the existence of two close 



1988), i.e., the existence of the 



bands in perovskites-like structures ([Sinha and Sinhal Il964l )). 



y situated filled p and d empty 



Eq. (|15|) has been used to calculate the electron-phonon contribution to the ferroelec- 
tric soft mode at the T point in PbTiC>3 and NaNbC>3. More precisely, we obtained the 
eigenvalues of the LDA-determined (q) matrix corresponding to the LDA-determined 
eigenvectors of the total dynamical matrix. The resulting eigenvalue iJ 1 corresponding to 
the soft TO mode is negative for any pressure and shown in Fig. EH and clearly behaves in a 
nonmonotonic way with pressure. Interestingly, this nonmonotonic behavior varies in a cor- 
related way with the soft-mode total frequency (see Figs. El El El) under pressure- which thus 
confirms that the electron-phonon coupling is the leading mechanism for the high-pressure 
ferroelectricity. 



IV. CONCLUDING REMARKS. 



In summary, our results reveal that, unlike commonly thought, ferroelectricity is not 
suppressed by high pressure in insulating perovskites. Instead, ferroelectricity is found to 
enhance as pressure increases above a critical value. Moreover, this unexpected high-pressure 
ferroelectricity is different in nature from conventional ferroelectricity because it is driven 
by an original electronic effect rather by long-range ionic interactions. Depending on the 
energetic difference between ferroelectric and paraelectric phases at atmospheric pressure, 
such phenomenon can lead to various overlooked/original effects, e.g. (i) the disappearance 
and then re-entrance of ferroelectricity under pressure; (ii) the occurrence of FE at high 
pressure in a nominally-paraelectric compound; and (iii) the existence of FE at any pressure 
(see Figure 6). We thus hope that this article, and its numerous details, will stimulate the 
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investigation of novel effects in "smart" materials under pressure. 
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FIG. 1 First-principles prediction of the pressure behavior of (a) ATI (see text) for all the considered 
phases (the phase corresponding to a minimal ATt at a given P is thus the most stable one for 
this pressu re), (b) the Cartesia n components (along the pseudo-cubic (001) axes) of the polar local 



soft-mode (jZhong et al 



1995) - which is directly proportional to the spontaneous polarization 



for the P4mm, I4cm, R3m and R3c phases, and (c) the rotational angle of the oxygen octahedra 
with respect to the pseudo-cubic (111) direction for the R3c and R3c rhombohedral phases and 
with respect to the pseudo-cubic (001) direction for the I4cm and I4/mcm tetragonal phases. 
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FIG. 2 The LDA-predicted pressure behavior of the c/a distortion at OK for the P4mm phase in 
PbTi0 3 . 
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FIG. 3 The pressure dependence of the elastic stiffness coefficients en (squares), cu (dots) , and 
C44 (triangles) of the cubic phase of PbTiOs, as predicted by LDA. 
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FIG. 4 Pressure evolution of the normalized Born effective charges for Pb, Ti, and On in the 
cubic phase of PbTiOs. The normalization is done with respect to the Born effective charges at 
zero pressure obtained with ABINIT for Pb, Ti, Oj_, and On in the cubic Pm3m phase - that are 
3.915, 7.14, -2.61 and -5.85, respectively. 
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FIG. 5 Pressure behavior of phonons in PbTiO"3 (solid symbols) and BaTi03 (open symbols). 
Panel (a) displays the square of the zone-center TO phonon frequency in the cubic phase within 
LDA at OK. Panel (b) shows the contributions Sflmi the Coulomb interactions (squares) within the 
rigid ions model and non-Coulomb interactions (triangles) to this square. 
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FIG. 6 LDA-predicted pressure behavior of the soft-mode eigenvalue and eigenvectors of the 
second-derivative matrix in cubic PbTiOs (panels (a) and (d)), BaTiC>3 (panels (b) and (e)) and 
BaZrO"3 (panels (c) and (f)). The GGA results are also shown for PbTiOs in panels (a) and (d), 
and suggest that PbTiOs may be paraelectric within a small pressure range. O3 is the oxygen 
atom located between two (slightly displaced) B-atoms along the [001] direction. 
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FIG. 7 LDA-predicted pressure behavior of the soft-mode eigenvalue and eigenvectors of the 
second-derivative matrix in cubic KNbO"3, NaNbO"3, SrTiO"3, CaTi03 and PbZr03. 
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FIG. 8 LDA-predicted pressure behavior of the soft-mode frequency square in cubic WO3. 
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FIG. 9 LDA-predicted pressure behavior of the soft-mode frequency square in the cubic phase of 
PbTiO-3, when excluding the 5d semicore states of Pb and 3s and 3p semicore states of Ti from 
the valence electrons (unlike the 6s and 6p states of Pb, the 3d and 4s states of Ti and the 2s and 
2p states of oxygen that are kept in the valence). 
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FIG. 10 LDA-predicted pressure behavior of the soft-mode eigenvalues of the second-derivative 
matrix in cubic bismuth aluminate BiAlC>3 and magnesium silicate MgSi03. 
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FIG. 11 Electronic-related properties of PbTiO"3 under pressure. Panels (a), (b) and (c) display the 
partial electronic density of occupied states in rlfie cubic phase at 3.6, 40 and 103 GPa, respectively, 
for the O 2s, O 2p and Ti 3d orbitals. The zero in energy is chosen at the top of the valence band, 
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FIG. 12 Integrated DOS of PbTi03 under pressure in the cubic state. Panel (a) displays the 
integrated PDOS for Ti 3d and O 2p and 2s states; panel (b) displays the normalized integrated 
PDOS for the same states; panel (c) displays the contributions to the normalized integrated PDOS 
coming from the low-lying O 2s and Ti 3d e g states (located between -22 and -15 eV) and from 
the topmost group of the valence O 2p and Ti 3d bands. 
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FIG. 13 Integrated PDOS of BaTiO-3 under pressure in the cubic state. The inset shows the 
variation of the integrated PDOS of Ti 3d e g and Ti 3d t2 g versus pressure. 
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FIG. 14 The bond degree parameter BD = h(r cp ) / p(r cp ) vs. \v(r cp )\/g(r cp ) ratio for the Ti-0 
(squares) and Pb-0 (triangles) bonds. The arrows indicate the direction of pressure increase. 
Vertical dashed lines indicate the regions of stability for the three classes of atomic interactions 
(see text). 
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FIG. 15 (a) Pressure evolution of the ionic Bader volumes for Pb, Ti and O and the unit cell 
volume \ 'cell in the cubic phase of PbTiO"3. (b) The pressure dependence of the normalized 
volumes. (The normalized volumes are equal t43unity at GPa.) 
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FIG. 16 LDA-predicted deviations from the Cauchy relation (5 = Cn — C44 — 2P) as a function 
of pressure for the cubic phase of PbTiO-3. 
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FIG. 17 The Ti3d-02s and Ti3<i-02p bonds polarity and covalency of a perovskite as a function 
of the interatomic Ti-0 distance d in the ideal cubic phase. (Note that the pressure increases as d 
decreases) . 
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FIG. 18 Change in bonding energy (solid lines) within the bond-orbital model of Harrison. This 
change is directly proportional to the soft-mode frequency square within this model. The contri- 
butions of the sd and pd bonds to that total change are also indicated via dashed lines, d is the 
Ti-0 distance in the ideal cubic phase. 
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FIG. 19 Normalized electron-phonon contribution uj^ to the ferroelectric soft mode at the T point 
in the cubic state of PbTi03 and NaNbC>3 . Normalization has been done such as to obtain values 
of -1 at GPa. 47 



